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The LHC Inverse Problem, Supersymmetry, and the ILC 
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We address the question whether the ILC can resolve the LHC Inverse Problem within the frame- 
work of the MSSM. We examine 242 points in the MSSM parameter space which were generated at 
random and were found to give indistinguishable signatures at the LHC. After a realistic simulation 
including full Standard Model backgrounds and a fast detector simulation, we find that roughly only 
one third of these scenarios lead to visible signatures of some kind with a significance > 5 at the 
ILC with ^/s = 500 GeV. Furthermore, we examine these points in parameter space pairwise and 
find that only one third of the pairs are distinguishable at the ILC at 5cr. 



The exploration of the Terascale begins in earnest next 
year with the start of operations at the LHC. Although 
the Standard Model (SM) does a respectable job at de- 
scribing present data, we know that it is incomplete for 
many reasons. One can argue persuasively that to ad- 
dress these issues new physics must appear at the Teras- 
cale. Once new physics is found, our next goal will be to 
determine the nature and detailed structure of the under- 
lying framework from which it arose. Can this be done 
uniquely using data from the LHC? The answer to this 
question has been recently quantified by Arkani-Hamed, 
Kane, Thaler and Wang (AKTW) Q, who demonstrated 
what has come to be known as the LHC Inverse Prob- 
lem. These authors studied a restricted and very spe- 
cific 'simple' form of new physics given by the Minimal 
Supersymmetric Standard Model (MSSM). By perform- 
ing a scan of MSSM parameter space, they showed that 
specific signatures observed in LHC experiments can not 
be uniquely mapped back to distinct parameter space 
points, hereafter referred to as models, within the MSSM. 



*This work is supported in part by funds provided by the U.S. De- 
partment of Energy (D.O.E.) under cooperative research agreement 
DE-FC02-94ER40818. 

■I-Work supported in part by the Department of Energy, Contract 
DE-AC02-76SF00515. 

ttResearch supported in part by the US Department of Energy 
under contract DE- AC02-06CH11357 

t Electronic address: pfberger@mit . edu 

5 Electronic address: ligaincr@slac. Stanford. edj 
^Electronic address: 



The purpose of the present Letter is to outline an inves- 
tigation of the capability of the International Linear Col- 
lider (ILC) to resolve this issue of degeneracy within the 
MSSM framework arising from the LHC Inverse Problem. 
Such an analysis, though of interest in its own right, pro- 
vides a unique opportunity to make a detailed study of 
the signals and backgrounds for hundreds of Supersym- 
metric (SUSY) models. From this we can re-examine our 
basic assumptions and prejudices about SUSY analyses 
at the ILC. 

An outline of our procedure is as follows^: (i) We ob- 
tained the set of model pairs from AKTW which led 
to identical signatures at the LHC. After some filter- 
ing, this consists of 242 models in 162 degenerate pairs 
with some models appearing in several pairs, {ii) We 
next generated signal (S) 'data' via Monte Carlo tech- 
niques for each of these models employing both PYTHIA 
[^ and CompHEP jj] and employing a realistic beam- 
spectrum [g specific to the ILC design. We assumed 
an ILC with a center of mass energy of -/i = 500 GeV 
and an integrated luminosity of 500 fb~ which was split 
equally between samples with 80% left- and right-handed 
electron beam polarization. Positron polarization was 
not included, (iii) We obtained two, statistically inde- 
pendent SM background sets (B) which were generated 
by T.Barklow [p| and employ full tree-level matrix ele- 
ments using WHIZ ARD/0 'Mega 0. These included all 
(> 1000) processes of the form 2 ^ 2, 2 ^ 4, and 2 ^ 6 
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^ For full details of our analysis see Ref. [E 
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TABLE I: Number of models, out of 242, which have a given 
final state kinematically accessible at ^/s — 500 GeV and 1 
TeV. Note that for the 500 GeV machine 96/242 models have 
only LSP or neutral pairs accessible while 61/242 models have 
no SUSY particles accessible. 



sions. The use of full matrix elements here is important 
as they generally produce larger backgrounds with longer 
tails in kinematic distributions than does, e.g., PYTHIA. 
(iv) We piped both signal and background through the 
java-based SiD B fast detector simulation, org.lcsim M, 
to take realistic detector effects into account. The pro- 
cedure outlined above corresponds to a realistic scenario 
of a first scan of the MSSM signature space at the ILC. 

At this stage, we first ask what are the numbers and 
types of SUSY particles that are kinematically accessible 
at the ILC in each of the 242 AKTW models? The an- 
swer is presented in Table |; here we see that many more 
sparticles are accessible at 1 TeV than at 500 GeV, which 
is a good argument for going to higher energies as soon 
as possible. 

The next issue to address is how many of the kinemat- 
ically accessible sparticles in a given model actually pro- 
duce an observable signal at the ILC. Certainly we can- 
not distinguish models which are invisible! To this end, 
we performed a number of detailed analyses to search 
for the selectrons, smuons, staus, sneutrinos, hghtest 
charginos and the lightest two neutralinos in each of the 
AKTW models, the details of which are given elsewhere 
PI . An important point here has to do with choosing the 
kinematic cuts to enhance the signal over background; 
many such sets of cuts have been developed in the liter- 
ature for ILC SUSY analyses ||l^; O. These are usually 



based on a few representative benchmark models, such as 
SPSla' [|2, 13|, which have very large SUSY production 
cross sections and many spartners kinematically accessi- 
ble. Since the AKTW models arise from a generalized 
parameter scan, these existing cuts are not always ade- 
quate for observing our SUSY signals, e.g., in some cases 
the existing cuts completely removed our signal, while in 
other cases they did not sufficiently reduce the SM back- 
ground. We thus developed a set of cuts that work more 
or less uniformly throughout the parameter space. This 
led in many cases to increased backgrounds compared 
with those found for the specialized points studied in the 
literature. 

In order to ascertain the significance, S, of the model 
signal for any given search analysis we employ the like- 
lihood ratio method, based on a Poisson statistics likeli- 
hood function, and use both of our statistically indepen- 
dent background sets as input. As is traditional, to claim 
that sparticle production is observable, we require that 
5>5. 




FIG. 1: Muon energy distribution: the number of events/2 
GeV bin (combined signal and background) after imposing 
selection cuts for several representative models, with right- 
handed electron beam polarization, assuming an integrated 
luminosity of 250 fb~^. For comparison, SPSla' is shown as 
the histogram with the largest number of events. The SM 
background is represented by the black histogram. 

An example of one of our analyses is shown in Fig. n. 
Here we display the single muon energy distribution from 
the production and decay of smuon pairs at the ILC; the 
specific process of relevance here is e+e^ — > )J-^ jj-^ — *■ 
^+/i^ -I- 2xi. Since here Xi is the lightest SUSY particle 
(LSP) and is both stable and neutral, it appears in the 
detector as missing energy and thus only the properties 
of the final state muons can be measured experimentally. 
Two things are apparent from this figure: (i) Although 
the backgrounds are relatively low due to judicious cuts, 
there is a wide range for the ratio S/B as the MSSM 
parameter space is scanned. Beam polarization plays an 



important role here, (ii) S/B for SPSla' is significantly 
larger than in any of the hundreds of models in our sam- 
ple; this result is generally found to hold in all of our 
other sparticle search analyses as well. 
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FIG. 2: Distribution of the chargino-LSP mass difference 
versus the chargino mass for 51 of the AKTW models with 
Am,^ < 6 GeV for the Xi states that are kinematically ac- 
cessible at ^/s — 500 GeV. The blue crosses represent models 
that are observable in our suite of analysis channels based on 
the Xi decay via off-shell W bosons. The green crosses corre- 
spond to models that are only visible in the radiative chargino 
production analysis channel, while the magenta ones represent 
models that yield observable signals in both the radiative and 
off-shell W channels. The black crosses are models that are 
visible in the stable chargino analysis. The red points are the 
4 models where the Xi state is not observable in any of our 
analysis channels, essentially due to phase space restrictions. 



The situation for other spartners, such as the lightest 
chargino Xi i is rnore complicated as the observable sig- 
nature strongly depends on the mass splitting with the 
lightest neutralino, Am = m-+ — m^o. If Atti is suffi- 

'^ ' Xl Al 

ciently large, one can observe the decay of pairs of x^'^ 
through on- or off-shell W^s into jets and/or leptons. For 
smaller values of Am, these decay products are soft, so we 
must tag the final state by requiring the associated pro- 
duction of a high-pT photon ||l^, ^ . When Am < m-,^ , 
Xi is effectively stable, decaying outside the detector, 
and we perform a search for long-lived charged particles 
|16| . By combining these various techniques, we can cover 
the entire MSSM parameter space fairly well as shown in 
Fig. g, where we find that only 4 out of the 53 models 
with accessible Xi are not observable at the ILC. 

Another potentially difficult case is stau pair produc- 
tion, e+e" — » f+f^ -^ T^T^ + 2x1 since the final state 
r's themselves decay and need to be identified in the de- 
tector. In our analysis, we focus on the hadronic decays 
of taus into pions, r -^ tti/,-; t — > pi^T- -^ ii tt^Vt] t — > 
Stti't-, the latter being a 3-prong jet, and also include the 
e, /i leptonic decays of the r as these make up a substan- 
tial branching fraction. As an alternative possibility, we 



allow leptonic tau decays into muons, but reject taus that 
decay into electrons. This significantly reduces contam- 
ination from photon-induced backgrounds, at the price 
of reducing the signal by roughly 30%. The reason for 
rejecting events with electrons is to reduce contamina- 
tion from gamma-induced Standard Model processes not 
involving taus. The following main background, which 
has large missing energy, and thus mimics stau pro- 
duction, is removed by the alternative tau ID method: 
6677 ~^ e(6)/i(/x) -I- E"''""', where the not reconstructed 
final states are indicated by the brackets ( ) , and a beam 
electron or positron is falsely identified as a tau decay 
product. Better forward muon ID capabilities would rem- 
edy this situation because the second muon would then 
be captured by the detector as well. 




FIG. 3: Reconstructed r energy distribution: the number 
of events/2 GeV bin (combined signal and background) af- 
ter imposing selection cuts for several representative models, 
with right-handed beam polarization, assuming an integrated 
luminosity of 250 fb^^. The bottom panel corresponds to re- 
moving the electron final state in r decays in order to avoid 
beam backgrounds as discussed in the text. The SM back- 
ground is represented by the black histograms. 

The results from both the standard and alternative f 
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TABLE 11: Number of models, at ^s = 500 GeV, which have 
a given final state particle visible above the SM background 
with a significance 5 > 5 divided by the number of models 
with the same particle kinematically accessible. 

analyses are shown in Fig. Q. While stau production is 
certainly observable above background in both cases, the 
alternative technique nearly completely eliminates the 
background. Augmenting the detector with muon identi- 
fication (ID) capabilities at lower angles could reduce the 
7-induced background without having to pay the price of 
introducing a restricted tau identification. Apparently a 
significant portion of the AKTW models have relatively 
low stau signal rates, and an improved muon ID capa- 
bility can be crucial if in fact this portion of the SUSY 
parameter space is realized in nature. Note that while 
21/28 models with kinematically accessible f's are vis- 
ible using the alternative ID, the number is reduced to 
only 13/28 if the standard method, which includes elec- 
tron final states, is employed. 

Table O summarizes the results of our individual model 
searches. Our analyses do well at finding charged spart- 
ners but not so well with neutral ones. We find that 
78/85 models with at least one charged spartner kine- 
matically accessible are visible, but only 17/96 models 
with only neutral spartners kinematically accessible are 
observed. This translates into 82 visible models out of 
161 models with kinematically accessible spartners, and 
82 models that are visible out of the total 242 AKTW 
models. 

Finally, now that we know which models lead to ob- 
servable signatures at the ILC, we can ask how well 
the ILC performs at distinguishing pairs of models that 
gave degenerate signals at the LHC. To this end we per- 
form a x^ comparison of the signal plus background 
histograms for each of our analyses for each of the 
model pairs. Specifically, we construct a x^ of the form 
X^iSl + B1,S2 + B2) where 5*1(2) are the pure signal 
samples for the two models being compared, and 51(2) 
are our two independent background samples. We find 
that we can distinguish, at the 5(3)cr confidence level, 
57(63)/72 pairs of models where at least one of the mod- 
els has a charged spartner kinematically accessible, but 



we fail completely when both models being compared 
have only neutral spartners accessible. Thus, out of all of 
the AKTW model pairs, our results show that 57(63)/162 
can be distinguished at the ILC 5(3)cr. 

From this analysis it is clear that the ILC with the 
SiD detector does a respectable job in observing charged 
sparticles that are kinematically accessible, and in distin- 
guishing models that contain such particles. The major 
weakness, beyond the restricted kinematic reach, is in the 
neutral spartner sector. This problem might be resolved 
by employing positron beam polarization [|l7[ as well as 
by implementing more sophisticated analyses. In order 
to observe more sparticles in all models, an upgrade to 
y/s = 1 TeV as soon as possible is desirable. 

We are indebted to Tim Barklow for extraordinary help 
with the Standard Model background and for many in- 
valuable discussions. 
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